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Abstract 

We review and extend some recent investigations of the effects of 
aperiodic interactions on the critical behavior of ferromagnetic g-state 
Potts models. By considering suitable diamond or necklace hierarchi- 
cal lattices, and assuming a distribution of interactions according to a 
class of two- letter substitution rules, the problem can be formulated in 
terms of recursion relations in parameter space. The analysis of stabil- 
ity of the fixed points leads to an exact criterion to gauge the relevance 
of geometric fluctuations. For irrelevant fluctuations, the critical be- 
havior remains unchanged with respect to the uniform systems. For 
relevant fluctuations, there appears a two-cycle of saddle-point char- 
acter in parameter space. A scaling analysis, supported by direct 
numerical thermodynamic calculations, shows the existence of novel 
critical universality classes associated with relevant geometric fluctu- 
ations. Also, we show that similar qualitative results are displayed by 
a simple model of two directed polymers on a diamond hierarchical 
structure with aperiodic bond interactions. 

PACS numbers: 64.60.Ak; 05.50.+q; 61.44.Br; 83.80.Rs 

1 Introduction 

In some recent publications [|T], ^, ^ we investigated the critical behavior 
of a g-state ferromagnetic Potts model on a class of hierarchical diamond lat- 
tices with aperiodic interactions. Taking advantage of the lattice structure, 
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and assuming a (layered) distribution of exchange interactions according to 
suitable two-letter substitution rules, we were able to write exact recursion 
relations in order to characterize the critical behavior under the influence 
of geometric fluctuations. For disordered ferromagnetic systems, the well- 
known Harris criterion predicts a change in critical behavior if the under- 
lying uniform system displays a positive speciflc heat critical exponent. An 
extension of this heuristic criterion has been proposed by Luck^] in order 
to gauge the relevance of geometric fluctuations associated with aperiodic- 
ity. We derived an exact expression of an analog of Luck's criterion on a 
hierarchical lattice. For weak fluctuations, there is a nontrivial symmetric 
fixed point in parameter space that leads to the same critical behavior of the 
uniform model. For strong enough fluctuations, this symmetric flxed point 
becomes fully unstable, and there appears a cycle-two in parameter space, 
that is shown to lead to a novel class of (aperiodic) critical behavior [g, ^. 
As the systems are quite simple, it has been feasible to perform a direct 
thermodynamic analysis of the free energy to check the predictions of scaling 
arguments. 

Now we review and extend these calculations. In Section II, we show 
that the same results are obtained for a g-state Potts model on a necklace 
hierarchical lattice. Although the critical temperature may change, critical 
exponents associated with either the symmetric fixed point or the cycle-two 
attractor do display universal features that do not depend on the particular 
choice of diamond or necklace hierarchical structures (and on particular val- 
ues of the interaction energies). In Section III, we consider a model of two 
directed polymers on a diamond lattice |P, with a (layered) distribution of 
aperiodic interactions, according to some two-letter substitution rules. The 
exact recursion relations in parameter space turn out to be very simple, but 
with a similar structure as in the case of the aperiodic Potts model. We 
then present some examples to show that, for relevant geometric fluctua- 
tions, a symmetric flxed point becomes fully unstable and there is again a 
novel cycle-two attractor of saddle-point character in parameter space. 

2 Potts model with aperiodic interactions 

The successive application of the period-doubling two-letter rule, 

A-^ AB, B ^ AA, (1) 
on an initial letter A, produces the sequence 

A^ AB ^ ABAA ABAAABAB .... (2) 
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At each stage of this construction, the numbers A^i and Nb-, of letters A and 
can be obtained from the recursion relations 

where the eigenvalues, Ai = 2 and A2 = —1, of the substitution matrix M 
govern most of the geometric properties of this sequence^. At order n, we 
can write 

N^A^ = \\1 + \\l and < = ^A^-iA^, (4) 

from which we have the asymptotic expressions A^]^"'' ~ A" and AA^^"'' ~ Ag , 
for large values of n, where AA^^"^ can be regarded as a fluctuation around 
an asymptotic value. As the total number of letters, at order n, is given by 
iV^"-) = A^j^"-* + A'^"\ we have AN^^ ~ (A^^"))'^, with the wandering exponent 

InlAal , . 

The ferromagnetic g-state Potts model is given by the Hamiltonian 

(id) 

where (Tj = 1,2, ...,q, at all sites of a lattice, Jij > 0, and the sum refers 
to nearest-neighbor pairs of sites (for g = 2, we regain the standard Ising 
model). We now consider a Potts model on a hierarchical necklace lattice, 
with 6 = 2 bonds and m branches, and assume that the couplings can take 
just two values, Ja and Jb, associated with a sequence of letters produced by 
the period-doubling substitution. In Fig. 1, we indicate some stages of this 
construction (for 6 = 2 and m = 3). Note that the period-doubling sequence 
is perfectly suitable for a hierarchical lattice with 6 = 2 bonds. Also, note 
that the choice of the same interactions along the branches turns out to 
mimic an aperiodic layered structure in the corresponding Bravais lattice. 

Now we decimate the internal degrees of freedom of the necklace structure 
and write exact recursion relations for the reduced coupling parameters. 
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Figure 1: Some generations of a hierarchical necklace lattice with 6 = 2 bonds 
and m = 3 branches. Letters A and B indicate the layered distribution of 
exchange interactions. 



where xa,b = exp {PJa,b) with j3 = l/ksT. 

In parameter space, besides the trivial fixed points, associated with zero 
and infinite temperature, there is always a non-trivial symmetric fixed point, 
associated with the critical behavior of the underlying uniform model (see 
the sketch in Fig. 2a). The linearization of the recursion relations about the 
symmetric fixed point x* leads to the matrix form 

where C (x*) is a structure factor and M is the transpose of the substitution 
matrix. The eigenvalues of this linear form are given by 

Ai = AiC(a;*) = 2C(a;*), (10) 

and 

= {x*) = -C {x*) . (11) 

It is easy to show that Ai > 1. Note that we can write Ai = 2^^, with 
yrp = D / {D — Ou), where D = In (2m) / In 2 is the fractal dimension of the 
lattice and a„ is the specific heat critical exponent of the underlying uniform 
model. Also, it is easy to show that IA2I < 1, for q < Qc, and IA2I > 1, for 
q > qc- In this last case, for q > qc, the symmetric fixed point becomes fully 
unstable, and the geometric fluctuations are relevant (see the sketch in Fig. 
2b). For b = m = 2, both for necklace and diamond structures, we have 
qc = 4: + 2^/2 = 6.828427..., which also corresponds to the critical number of 
states for the relevance of disorder in a ferromagnetic Potts model on a fully 
disordered diamond lattice [Q. 
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(a) (b) 
Figure 2: Flow diagrams in the neighborhood of the symmetric fixed point. 



For necklace and diamond lattices with b bonds and m branches, and a 
suitable two-letter substitution sequence with period Ai = b, we can always 
write 



Ai =2C (x*) = 6^/(2-""), 



:i2) 



where D = In (m6) /In 6 is the fractal dimension of the lattice, is the 
specific heat critical exponent of the underlying uniform model, and 



A, 



A2C (x*) = X'^C (x*) = (x*) 



where the wandering exponent uj comes from Eq. 



We then have 



A, 



(13) 



(14) 



that leads to an exact version of Luck's criterion |P for the relevance of geo- 
metric fiuctuations. 



u; > 1 - 



D 



a,. 



(15) 



For the g-state Potts model on both diamond and necklace lattices with 
6 = 2 bonds and m = 2 branches, this inequality is equivalent to g > gc = 
4 + 2^2, as we have pointed in the last paragraph. In this special case, it 
also corresponds to > 0, which is in agreement with the usual form of 
Harris criterion for the relevance of disorder 0]. It should be pointed out 
that, according to the usual expectations, given q, b and m, the eigenvalues 
Ai and A2 assume the same values for both diamond and necklace structures. 

We now turn to the investigation of the critical behavior when the sym- 
metric fixed point is fully unstable. Numerically, we can see that, instead of a 
hypothetical aperiodic fixed point, there appears a two-cycle, with a saddle- 
point character, in parameter space. We then invoke scaling arguments to 
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write 



f{x)=g (x) + --rf (x") , 



(16) 



where / (x) is the reduced free energy per bond, g (x) is a regular function, 
and x" is a second iterate of the recursion relations. Note that the factor 
has to be included because we need two iterates to go back to the vicinity of 
the initial point in parameter space. This equation leads to the asymptotic 
solution |[T0[] 



where x* is one of the points of the two-cycle, Acyc is the largest eigenvalue 
of the linearization of the second iterate of the recursion relations about any 
one of the points of the cycle, P [z) is an arbitrary function of period 1, and 
the specific heat critical exponent is given by 



Let us show some calculations for the g-state Potts model on a neck- 
lace lattice, with 6 = 2 bonds, m branches, and a distribution of exchange 
interactions according to the period-doubling two-letter sequence: (i) for 
q = 7 and m = 2, the two-cycle is located at {x*j^,x*^) = (2.299..., 2.764...) 
and (2.587..., 2.179...), with eigenvalues of the second iterate given by Ai = 
3.993... and Aa = 0.985.... From Eq. (|T8|), we have a = -0.00226..., to 
be compared with the positive value a„ = 0.010..., for the uniform sys- 
tem; (ii) for g = 25 and m = 2, the two-cycle is located at (x^,x^) = 
(2.634..., 15.308...) and (6.246..., 1.957...), with eigenvalues of the second it- 
erate given by Ai = 4.243... and A2 = 0.343.... From Eq. (0), we have 
a = 0.08174..., to be compared with = 0.40456..., for the corresponding 
uniform system. For other values of m, we find qualitatively similar results. 
The specific heat exponent is definitely depressed as compared to the values 
for the underlying uniform system. Also, given q, b, and m, the exponents 
assume the same values for diamond and necklace structures. 

To check the validity of the scaling arguments, and the role of the two- 
cycle as the responsible for the new critical behavior, we have performed 
direct numerical analyses of the singularity of the thermodynamic free-energy. 
For the g-state Potts model on the 6 = 2 necklace hierarchical lattice, and 
with the period-doubling rule, we can write the free energy per bond as a 




(17) 





6 



series expansion, 



f{xA,XB] 



n=0 



3m 
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in) 
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+ q-2 



H In 

6m 



2 X 



In) 



+ q 



(19) 



All of our numerical checks fully confirm the scaling results. For example, for 
m = 2 and q = 100, with Ja/ Js = 5, a numerical analysis of the specific heat 
divergence at the critical temperature leads to a = 0.27 ± 0.03, to be com- 
pared with the scaling prediction from the two-cycle, a = 0.27204... (which 
is definitely different from the value for the uniform system, = 0.64846...). 
Within the error estimates, we always obtain the same numerical results for 
all values of the ratio Ja/Jb- Again, given g, 6, and m, and for all val- 
ues of the ratio JaIJb-, the exponents assume the same values for diamond 
and necklace structures. From the numerical calculations, and particularly 
for bigger values of g, we have been able to detect a log-periodic oscillatory 
behavior of the thermodynamic functions. The period of oscillation is com- 
patible with Eq. (17), which has also been confirmed in the independent 



calculations of Andrade[|TT|. 



3 Directed polymers with aperiodic interac- 
tions 

Along the lines of a paper by Mukherji and Bhattacharjee[P], we now consider 
a model of two interacting directed polymers on a diamond lattice (see Fig. 
3). The polymers start at one end of the lattice and meet at the other end. 
There is an attractive interaction if a bond of the lattice is shared by two 
polymers. 

Consider a diamond lattice (6 = 2 bonds and m branches), with a lay- 
ered distribution of interactions, according to the two-letter period-doubling 
sequence given by Eq. (|1|). It is easy to write the recursion relations 

ffi — 1 

v'a = -VaVb + , (20) 

m m 

and 

, 1 ^ m — 1 , , 

yB = -y\ + , 21 

m m 
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Figure 3: Diamond hierarchical lattice with b = 3 bonds and m = 2 branches. 



where yA,B = exp {Pva,b), and va,b > is the interaction energy at bonds 
of types A and B, respectively. For i/a = ys = y, recover the recursion 
relation for the uniform model, 

, 1 9 m — 1 , , 

y' = -y' + • 22 

m m 

Besides the trivial fixed points, y* = 1 and cxo, associated with zero and 
infinite temperatures, there is a nontrivial fixed point, y* = m — 1, which 
becomes physically acceptable, and is associated with a binding-unbinding 
transition, for m > 2 (there is no phase transition on the simple diamond 
lattice with m = 2 branches). 

The recursion relations ( |2D| ) and ( |?T| ) are so simple, that it easy to see 
that, for m > 2, there is no physically acceptable nontrivial fixed points 
except the symmetric fixed point y*A = y*B = y* The linearization of the 
recursion relations in the neighborhood of this symmetric fixed point leads 
to matrix form 



Ay'^\^rfl l\(AyA 
Ay's J m\2 J \ Ays 



(23) 



with eigenvalues Ai = 2y*/m = 2 (m — 1) /m, and A2 = —y*/m = —{m — 
l)/m. Therefore, if m > 2, we have Ai > 1 and IA2I < 1, which shows that 
the geometric fluctuations are completely irrelevant in this case. 

We now turn to a more interesting case. Consider a diamond lattice with 
b = 3 bonds and m branches (in Fig. 3, we sketch a diamond lattice with 
b = 3 bonds and m = 2 branches). Suppose that the (layered) interactions 
are chosen according to the period-3 two-letter sequence, A — > ABB and 
B AAA. The substitution matrix is characterized by the eigenvalues 
Ai = 6 = 3 and A2 = —2, with the wandering exponent u = In 2/ In 3 = 
0.630092.... The new recursion relations are given by 

, 1 9 m — 1 , , 

yA = -VAyl + , 24 

m m 
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and 



, 1 o m — 1 

m m 

There is just a single nontrivial fixed point, for m > 3, at a symmetric 
location, 

1 1 



y*A^y%-y*^-^ + ^V^n^. (26) 

The linearization in the neighborhood of this fixed point leads to the matrix 
form 



^y's J m \S J \Ay, 
with the eigenvalues 



(27) 



y*2 3 

Ai = 3— = \2m-l- V4m-3| , (28) 

m 2m 



and 



*2 1 

A2 = -2— = \2m-l- V4m-3| . (29) 

mm 

For 3 < m < 3 + -\/5, it is easy to show that Ai > 1, and IA2I < 1. As in the 
case of the simple diamond lattice with 6 = 2 bonds, geometric fluctuations 
are irrelevant and the critical behavior is identical to the uniform case. How- 
ever, for m > 3+\/5 = 5.236068..., wc have IA2I > 1, and the symmetric fixed 
point becomes fully unstable. For example, for m = 5, we have y*A = y*B — 
y* = 1.561552..., with eigenvalues Ai = 2.140568... and A2 = 0.951363... < 1. 
For m = 6, however, we have y*A = y*B = V* = 1.791287..., with eigen- 
values Ai = 2.573958... and A2 = 1.143981... > 1. As in the case of the 
Potts model, there is then a cycle-two, of hyperbolic character, in parameter 
space. It is easy to locate this cycle at {y*A-,y*B) = (1.419001..., 2.267305...) 
and (2.049103..., 1.309541...), with eigenvalues of the second iterate given by 
Ai = 2.624300... > 1 and A2 = 0.772598... < 1. 



4 Conclusions 

In conclusion, we have taken advantage of suitable diamond and necklace 
hierarchical lattices, and used a distribution of (layered) interactions accord- 
ing to a class of two-letter substitution rules, to perform an exact analysis 



9 



of the effects of geometric fluctuations on the critical behavior of ferromag- 
netic g-state Potts models. We derived an exact expression for an analog 
of Luck's criterion on a hierarchical lattice. Also, due to the simplicity of 
the formulation, we have been able to show that the novel critical behav- 
ior under the influence of relevant geometric fluctuations is governed by a 
two-cycle attractor in parameter space. This new universality class can be 
described by scaling arguments, that are fully supported by direct thermo- 
dynamic calculations for the singularity of the free energy (as usual, critical 
exponents do not depend on details of the systems, as the relative strength of 
the interactions). It is not difficult to extend these calculations for an Ising 
model in the presence of direct and staggered fields[|l^, to analyze the effects 
of aperiodicity on the tricritical behavior in a mixed-spin Ising system[| 



and to consider some noncommutative sequences [|T^] 

We have also reported some results for the critical behavior of a model 
of two directed polymers on a diamond lattice with a (layered) distribution 
of aperiodic interactions (according to suitable two-letter substitution rules). 
The exact recursion relations in parameter space turn out to be very simple, 
but with a similar structure as in the case of the aperiodic Potts model. We 
have presented some examples to show that, for relevant geometric fluctu- 
ations, a symmetric fixed point in parameter space becomes fully unstable 
and there is again a novel cycle-two attractor, of saddle-point character, that 
governs a new class of (aperiodic) critical behavior. 

This work has been supported by the Brazilian agencies FAPESP and 
CNPq. 
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